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Abstract 


Campos-Díaz, K. E., £ Alvarez-Cruz, J. L. (May-June, 2017). 
A proposal of a hydrodynamic model to low Reynolds 
numbers in a liquid-solid inverse fluidized bed reactor. 
Water Technology and Sciences (in Spanish), 8(3), 143-150. 


This paper examines the fundamental characteristics of an 
inverse liquid- solid inverse fluidized bed reactor. In this 
reactor, many experiments were performed using three sizes 
of polypropylene spherical solids, lighter than water, with 
different diameters (4.00, 4.16 and 4.18 mm) and of densities 
ranging from 808 to 867 kg/m, these were fluidized by 
continuous liquid phase flow. Various bed heights were 
obtained at different flowrates from 0.95 to 9.5 (L min”). 
A predictive model is developed for the hydrodynamic 
behavior, such as height and solid holdups (bed porosity), 
this model was derived from the balance of forces acting 
on a single particle and the set of particles found in the 
fluidized bed; uses process parameters, such as reactor 
dimensions, particle properties, and liquid flowrates as 
input variables, this model also includes the inertial drag 
coefficient which allow extended a non-spherical particles. 
This model has several advantages when compared with 
previously reported. Among them, advantages such as 
standard deviation values <0.9% between experimental and 
calculated bed porosity, the bed porosity can be expressed 
in terms of the Reynolds and Archimedes numbers and this 
model estimate bed porosity in a range of Reynolds number 
from 5.5 to 200 which has not been studied in an inverse 
fluidized bed reactors. 
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Resumen 


Campos-Díaz, K. E., € Alvarez-Cruz, J. L. (mayo-junio, 2017). 
Propuesta de un modelo hidrodinámico para números de Reynolds 
pequeños en un reactor de lecho fluidizado inverso sólido-líquido. 
Tecnología y Ciencias del Agua, 8(3), 143-150. 


El presente estudio tuvo como finalidad estudiar el comportamiento 
hidrodinámico de un reactor de lecho fluidizado sólido-líquido 
operado a flujo inverso. En este reactor se realizaron varios 
experimentos usando partículas esféricas sólidas de polipropileno 
más ligeras que el agua de tres diferentes tamaños, con distintos 
diámetros (4.00, 4.16 y 4.18 mm) y densidades que van de 808 a 
867 kg/m, que fueron fluidizadas por un flujo continuo de líquido. 
La altura del lecho fue determinada a diferentes velocidades de 
flujo de 0.95 a 9.5 (1 min”). Se desarrolló un modelo para predecir 
el comportamiento hidrodinámico, como la altura y el vacío entre 
los sólidos (porosidad del lecho); este modelo se derivó de la relación 
de fuerzas que actúan sobre una sola partícula y el conjunto 
de partículas que se encuentran en el lecho fluidizado; utiliza 
parámetros del proceso, como dimensiones del reactor, propiedades 
de las partículas y velocidades del líquido como variables de entrada; 
este modelo también incluye el coeficiente de arrastre inercial, el 
cual permite aplicarlo para partículas no esféricas. Este modelo tiene 
varias ventajas comparado con los reportados por otros autores, entre 
ellas, valores de desviación estándar < 0.9% entre la porosidad del 
lecho experimental y la calculada; la porosidad del lecho puede ser 
expresada en términos de los números de Reynolds y Arquímedes, 
y este modelo estima la porosidad del lecho en un rango de número 
de Reynolds a partir de 5.5 a 200, el cual no ha sido estudiado en un 
reactor de lecho fluidizado inverso. 


Palabras clave: reactor de lecho fluidizado inverso, porosidad del 
lecho, porosidad, número de Reynolds, número de Arquímedes, 
hidrodinámica, vacío entre sólidos, coeficientes de arrastre inercial. 
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Introduction 


Fluidization is an operation in which solid 
particles enters in contact with a fluid acquir- 
ing fluid characteristics (Mukherjee, Mishra, éz 
Ran, 2009) and is known for its wide use since 
the early 50's. In conventional fluidization, 
the solid particles have a higher density than 
the fluid therefore the bed solids can be fluid- 
ized by an upwards flow, in the case of a bed 
of particles having a density smaller than the 
fluid (usually liquid), the bed is fluidized by a 
downwards flow of the liquid and it is usually 
called inverse fluidized bed (Bimal, Uma, éz Su- 
dip, 2010). Among the advantages presented by 
the latest are high mass transfer rates, minimum 
carry over of coated microorganisms due to less 
solid attrition, efficient control of biofilm thick- 
ness and ease of refluldization in case of power 
failure (Garcia, Buffiére, £ Elmaleh, € Moletta, 
1998; Bimal et al., 2010). These significant ad- 
vantages found many applications of inverse 
fluidized beds in biochemical processes like fer- 
rous iron oxidation and aerobic and anaerobic 
biological wastewater treatment like treatment 
of wine distillery wastewater. Previous studies 
on liquid-solid inverse fluidized bed reactor 
(IFBR) dealt with different hydrodynamic 
characteristics such as pressure drop, minimum 
fluidization velocity and bed void fraction e.g. 
(Fan, Muroyama, €: Chern, 1982; Karamanev 
£ Nivkolov, 1992; Renganathan € Krishnaiah, 
2007). 

The bed void fraction is one of the important 
design parameters that determine the height of 
fluidized bed during steady state operation. Un- 
der normal operating conditions, fluidized bed 
reactors operate under steady state. However, 
unsteady state operation is encountered during 
start-up and shut-down of the reactor and due to 
fluctuations in liquid flow rates. Such conditions 
may prevail in a wastewater effluent treatment 
plant where an IFBR is mostly used. Therefore, 
information on the unsteady state bed expan- 
sion is required for better understanding of 
the fluid dynamics under transient conditions, 
which will lead to better control of the reactor 
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and for the design of IFBR an industrial scale 
(Renganathan é Krishnaiah, 2007). 

In many processes there are particles with 
diameter different in the particles of the bed 
greatly modifies its porosity due to the particle's 
size and density variation. This modification in 
the bed porosity affects the proper operation 
and control of the fluid bed reactor. Therefore, 
itis important to have a reliable model in order 
to estimate bed porosity for design and scale-up 
of fluid bed reactors. 

Numerous models have been proposed to 
predict the bed porosity in fluidized beds with 
spherical and non-spherical particles, but only it 
includes a range of Reynolds number from 200 
to 6,000. The major and most used until today 
by several authors (Yang $ Renken, 2003; Akgi- 
ray éz Soyer, 2006; Renganathan éz Krishnaiah, 
2007; Fuentes, Scena, Aguirre, €: Mussati, 2008; 
Soyer $: Akgiray, 2009) are: Richardson and 
Zaki (1954), Wen and Yu (1966), Ramamurthy 
and Subbaraju (1973), Riba and Couderc (1977), 
Fan et al. (1982), and Setiadi (1995) models. 

Mathematical models for prediction of tra- 
ditional fluidized bed porosity can be used in 
inverse fluidization, nevertheless, some authors 
disagrees with this practice (Wen € Yu, 1966; 
Ramamurthy € Subbaraju, 1973; Hyun, 2001). 

Therefore the general aim of this work is 
to study the hydrodynamic of a liquid-solid 
inverse fluidized bed reactor with different 
polypropylene spherical particles diameters 
and different densities to propose a mathematic 
model to estimate bed porosity in a range of 
Reynolds number from 5.5 to 200, which has 
not been studied in this reactor and this hydro- 
dynamic characteristics are very important to 
low fluidization rates in biotechnology process. 


Materials and methods 

Liquid and solid phase 

Water at 20 *“C was used as fluid work in an 
inverse fluidization reactor and polypropylene 


spherical particles of average diameter 4.00, 4.16 
and 4.18 mm and density of 808, 825 and 867 
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kg/m? were used as solid phase. The sphericity 
is an important measure that describes how 
round a particle is defined as the ratio of the 
surface area of a sphere with the same volume 
as the given particle to the surface area of the 
particle itself, therefore the solid particles in this 
work have a sphericity of 1. 


Installation 


A schematic diagram of the fluidized bed reac- 
tor used in this work is shown in Figure 1. As 
shown, the experimental setup included a 0.04 
m of internal diameter (ID) and 0.50 m high 
cylindrical glass column. The work fluid was 
stored in a tank, and pumped into the column 
by a 1/8 HP little giant centrifugal pump, model 
583002. The flow rate in the reactor was mea- 
sured with a rotameter and fixed in the range 
between 0.95-9.5 (l min”). 


Start-up 
The inverse fluidization with each particles lot 


was done of the following way: The inverse 
fluidized bed reactor was filled with water at 


20 *C, the water was used as a liquid phase, the 
polypropylene particles were introduced for the 
upper part of the reactor and the initial height 
bed was measured with the calibrated scale in 
centimeter. The flow rate was increased 6 times 
more (0.63, 0.88, 1.00, 1.13, 1.26 and 1.38 1 min”) 
and the heights bed was measured. Figure 2 
shows the inverse fluidization reactor at differ- 
ent flowrates. 


Particle density 


Particle density was determined by Standard 
test methods for specific gravity of solids by wa- 
ter pycnometer (ASTM D854-14) to determining 
their volume and mass, this method was repeat 
three times The volume was measured using a 
pycnometer filled with light oil (density 750 kg/ 
m') and the mass was measured in an analytical 
balance Sartorius +1x107 kg. The particle den- 
sity was calculated by dividing the mass by the 
volume. 

The density media was determined by pyc- 
nometer (ASTM D854-14). The viscosity media 
was measured with a Malvern SV-10 vibratory 
viscosimeter. 


e 


wWJ GV 


Specifications 


BP = Bed of solid particles  R= Rotameter 
C = Column 

D = Distributor 

H = Homogenizer 
LI = Level indicator 
FT = Feed tank 

cp S = Scale 


WJ = Water jacket 
GV = Globe valve 
CV = Gate valve 
TCI = Temperature controller indicator 
CP = Centrifugal pump 
.... = Process water 
o - Inverse fluidization 


Figure 1. Schematic diagram of the inverse fluidized bed reactor. 
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Figure 2. Inverse fluidization reactor at different flowrates. 


Results and discussion 
Hydrodynamic study 


The different equations that relate the bed po- 
rosity with the fluid superficial rate can be clas- 
sified into three main groups Fan et al. (1982). 
The first group is based on a relationship that 
dependence between velocity and initial veloc- 
ity (U/U,) and porosity (e). The Richarson and 
Zaki (1954) equation is the most popular of this 
group. The second group of equations is based 
on the drag for multi-particulates systems. 
This is usually given as a function of Reynolds 
number (Re) and Archimedes number (Ar). 
Ramamurthy and Subbaraju (1973), and Riba 
and Couderc (1977) are typical for this group. 
The third group is based on the ratio ob- 
tained from the balance of forces acting on a 
particle isolated (gravity, floating, drag) and 
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assume that applies to all particles in the bed as 
the equation of Wen and Yu (1966) for spherical 
and non-spherical particles using a shape factor 
(yw) and Reynolds numbers are in the range 200 
to 6000. Table 1 shows the function of porosity 
shown Wen and Yu (1966), Riba and Couderc 
(1977), and Rammamurthy and Subaraju (1973) 
which were chosen because they have a lower 
standard deviation in the case of spherical par- 
ticles of uniform size. 


Predictive model proposed and comparison 
with other models 


The comparison of the results was made by plot- 
ting the porosity function obtained from values 
proposed by Wen and Yu (1966), Ramamurthy 
and Subbaraju (1973), Riba and Courdec (1977) 
and this work, versus experimental bed porosity 
data. 
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Table 1. Porosity function model by different authors. 


Author 


Porosity function 


fe) 
Wen and Yu (1966) Hwy = pa A 
WXY 18Re+2.7Re**% 
Riba and Couderc (1977) Mie = came 


1.9(18Re+2.7 Re!) 


Ramamurthy and Subbaraju (1973) 


1 


Ar F 


Hrs = laa 


Ar Archimedes number [- dSp,(p, -P5)/ 1] 


d, Particle diameter 
f Function 

fe) Porosity function 
2 gravity 


Ga Galileo number [- d;Spi/ wi] 
Mv Density ratio [= P,-P, /p,] 

u Liquid viscosity 

Re Reynolds number [= d,vp, /u ,] 
P Liquid density 

p, Particle density 


The model proposed in this paper was 
designed using the inertial drag coefficient (C) 
Becker (1959) for Re < 200. The data used in the 
model proposed in this work were Wilhelm 
and Kwauk (1948) and the experimental data 
obtained in this work. The model to calculate 
the porosity of spherical particles and the 
application intervals as shown in table 2. 

Ramamurty and Subaraju (1973) porosity 
function fits the experimental data tendency, but 
do not represent the experimental bed porosity. 
The plot (figure 3) shows an experimental un- 
derestimation for experimental data between 0.5 
and 0.73 of bed porosity, this deviation can be 
attributed to the principal assumption by these 
authors that the solid particles are represented 
as an ensemble moving about points considered 
as nodes of an imaginary lattice through whose 
free volume the fluid is flowing. Itis important 
to notice that this model was developed for an- 
nular fluid. 

Riba and Courdec (1977) model fits the 
experimental data tendency but the plot shows 
a slight experimental underestimation for ex- 
perimental data between 0.55 and 0.63 of bed 


porosity. 


Wen and Yu (1966) model fits has a good 
correlation of bed porosity data; the standard 
deviation was 1.2%. 

The mathematical model proposed in this 
work fits well the trend of the experimental bed 
porosity data for spherical particles which was 
well represented by this function in the inter- 
val of porosity from 0.47 to 0.73; the standard 
deviation was 0.9%. 

The new model has several advantages 
because Wen and Yu (1966), Riba and Cour- 
derc (1977), and Ramammurthy and Subaraju 
(1973), restrict their model to traditional fluidi- 
zation (upward flow) and spherical particles. 
Furthermore the new model is considerably 
simple in form, it can be used for spherical 
particles, but includes a term of inertial drag 
coefficient proposed by Becker (1959) which 
includes a shape factor to non-spherical 
particles. 

The figure 3 shows the graph of the func- 
tion of porosity that results from using the cor- 
relation Wen and Yu, Riba and Couderc (1977), 
Ramamurthy and Subbaraju (1973), and the 
model proposed in this work. 
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Table 2. Model proposed in this work to calculate the porosity of spherical particles. 


Model proposed in this paper 


Ar 


15 E 
0.75(24Re +C, Re?) 


Application intervals 


5.5 < Re < 200 
0.47 < € <0.73 
Spheres 


Ar Archimedes number [= disp, (p, -ps)/ ml 


Es Inertial coefficient 
d, Solid diameter 

E Porosity 

E Gravity 


. Liquid viscosity 

Re Solid Reynolds number [= d.vp, /u h] 
P. Liquid density 

P, Solid density 


fe) 


0.4 0.45 0.5 0.55 


——— Riba-Couderc equation 
=== This work equation 
= - — Wen-Yu equation 
x Experimental Wen-Yu 
+ Ramamurthy-Subbaraju equation 
Experimental Ramamurthy-Subbaraju 


Experimental Riba-Couderc 


Experimental this work 


0.6 0.65 0.7 0.75 0.8 


e (porosity) 


Figure 3. Model proposed to calculate the porosity of spherical particles. 


The table 3 shows a comparison of the stan- 
dard deviation of the equation proposed in this 
work and the equations published in literature. 


Conclusions 


The experiments were conducted in an inverse 
fluidized bed reactor at different flowrates. The 
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results show that the approach to the balance 
of forces acting on a particle, namely gravity, 
the buoyancy force and the drag applied to a 
normal flow, i.e., from the bottom up, is valid 
for case of a reverse flow (from top to bottom). 
The equation Wen and Yu (1966) and the equa- 
tion proposed in this work are showing a lower 
standard deviation to the experimental data. 
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Tabla 3. Standard desviation. 


Equation Standard desviation (%) 
Wen and Yu (1966) 1.2 
Riba and Courderc (1977) 2.69 
Ramammurthy and Subaraju ( 1973) 3.46 
Model proposed in this work (2016) 0.9 


The equation proposed in this paper can 
make a better approach by the balance of forces 
with a larger number of data in the literature. 
It is important to identify the study variables, 
also, to a dimensional analysis, since grouping 
variables in dimensionless numbers such as the 
Reynolds and Archimedes allow the handling 
of mathematical models takes just a simple 
manner, also, dimensionless numbers handling 
facilitates the design and scale reactors. 
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